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ABSTRACT 

We report the first observation of a transient relativistic jet from the canonical black 
hole candidate, Cygnus X-l, obtained with the Multi-Element Radio-Linked Inter- 
ferometer Network (MERLIN). The jet was observed in only one of six epochs of 
MERLIN imaging of the source during a phase of repeated X-ray spectral transitions 
in 2004 Jan-Feb, and this epoch corresponded to the softest 1.5-12 keV X-ray spec- 
trum. With only a single epoch revealing the jet, we cannot formally constrain its 
velocity. Nevertheless, several lines of reasoning suggest that the jet was probably 
launched 0.5-4.0 days before this brightening, corresponding to projected velocities of 
0.2c < v app < 1.6c, and an intrinsic velocity of > 0.3c. We also report the occurrence 
of a major radio flare from Cyg X-l, reaching a flux density of ~ 120 mJy at 15 GHz, 
and yet not associated with any resolvable radio emission, despite a concerted effort 
with MERLIN. We discuss the resolved jet in terms of the recently proposed 'unified 
model' for the disc-jet coupling in black hole X-ray binaries, and tentatively identify 
the 'jet line' for Cyg X-l. The source is consistent with the model in the sense that 
a steady jet appears to persist initially when the X-ray spectrum starts softening, 
and that once the spectral softening is complete the core radio emission is suppressed 
and transient ejecta / shock observed. However, there are some anomalies, and Cyg 
X-l clearly does not behave like a normal black hole transient in progressing to the 
canonical soft / thermal state once the ejection event has happened. 
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1 INTRODUCTION 

Relativistic jets (e.g. Hughes 1991; Burgarella, Livio & 
O'Dea 1993; Livio 1996) are a common feature of accretion 
onto compact objects, most notably black holes, on all scales. 
They are important channels for the removal of gravitational 
potential and/or spin energy, as well as angular momentum, 
and yet the details of their formation remain largely unclear. 
Recent studies of stellar-mass (3M© < Mbh < 15M©) black 
holes in binary systems have led to some empirical under- 
standing of the coupling between accretion 'states' and the 
jet formation process (Fender, Belloni & Gallo 2004 and ref- 
erences therein). Quantitative scalings between these low- 
mass black holes and the supermassive black holes in dis- 
tant galaxies (e.g. Merloni, Heinz & di Matteo 2004; Falcke, 
Kording & Markoff 2005) have finally given firm grounding 
to the hope that by studying X-ray binaries we may further 
our understanding of active galactic nuclei. 



Cygnus X-l (HDE 226868, V1357 Cygni) comprises a 
supergiant secondary (spectral type 09.7 lab) with mass 
between 20-33 M© , together with a compact primary which 
is a strong black-hole candidate with mass between 7-16 Mq 
(e.g. Gies & Bolton 1986a) . The system spends most of its 
time in a bright but hard X-ray state, with a bolometric X- 
ray luminosity around 2% of Eddington (e.g. di Salvo et al. 
2001). 

High angular resolution Very Long Baseline Array 
(VLBA) + phased Very Large Array (VLA) observations 
from 1998 Aug show a clearly resolved and persistent jet over 
three closely spaced epochs in the low/hard state (Stirling 
et al. 2001). This was the first direct imaging of steady ra- 
dio jet from a source in the hard X-ray spectral state. More 
recently, Gallo et al. (2005) have discovered a large-scale 
(^5 arcmin) radio and optical lobe approximately orien- 
tated with the mas-scale jet. Assuming the lobe is the result 
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Figure 1. Radio and X-ray monitoring of Cygnus X-l, from the Ryle Telescope and RXTE ASM. All data points arc daily averages. 
The top panel shows the radio and X-ray fluxes, and the lower panel the ratio of counts in the (5—12) / (1.5—3) keV bands; higher values 
indicate 'harder' X-ray spectra. Clear periods of radio:X-ray correlations (hard X-ray state) and anticorrelations (soft[er] X-ray state) 
are evident. The epochs of MERLIN observations are indicated by arrows at the top of the figure. Only the first MERLIN run, when the 
core radio source was very faint, revealed a jet. 



of the prolonged action of the jet on the ambient medium, 
Gallo et al. (2005) were able to demonstrate that the jet 
must be carrying a time-averaged power comparable to the 
present-day X-ray luminosity. 

2 OBSERVATIONS 

2.1 Ryle Telescope and RXTE ASM monitoring 

The Ryle Telescope (RT) is used from time to time for mon- 
itoring variable sources, including Cygnus X-l, at 15 GHz. 
The basic parameters of the RT and its monitoring pro- 
gram are described in Pooley & Fender (1997). The approx- 
imately daily monitoring of Cygnus X-l with the RT is com- 
plemented by data from the Rossi X-ray Timing Explorer 
(RXTE) all sky monitor (ASM). In addition to tracking the 
evolution of the total X-ray flux of the source, the RXTE 
ASM data also allow us to measure hardness ratios, indica- 
tive of the X-ray 'state' of the source (e.g. Fender, Belloni &: 
Gallo 2004; Remillard & McClintock 2005; Homan & Bel- 
loni 2005; Belloni 2005). Daily-averaged data from the RT 
and RXTE ASM, centred on the period of our MERLIN ob- 



servations (see below) are presented in Fig 1; see also Figs 
3-6. 



2.2 MERLIN 

MERLIN was initially triggered to observe Cyg X-l by a 
change in the X-ray behaviour of Cyg X-l to a generally 
brighter and more variable state, commencing around 2003 
December (see Fig 1). This trigger resulted in observations 
at epochs Ml and M2 (see Table 1). A major radio flare 
detected by the RT on 2004 Feb 20 resulted in a second 
MERLIN trigger and observations at epochs M3-6. We ob- 
served while the array was undergoing engineering work and 
so the specific telescopes observing changed between epochs. 
Altogether we obtained varying amounts of data over 8 sep- 
arate days. All of the MERLIN observations were taken at 
4.994 GHz and with 16 MHz bandwidth in each of the LL, 
RR, LR and RL polarization products. 

The epochs of the MERLIN observations are marked 
on Fig 1. From 2004 Feb 20 - 2004 Feb 22 Cambridge was 
involved sporadically (about 50, 10 and 30 per cent respec- 
tively of each days' observing time) and Knockin had a warm 
receiver cryostat throughout, it was removed from the array 
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Figure 2. A resolved transient jet from Cygnus X-l. This image 
corresponds to epoch Ml when the core was very faint, and is 
the only one of six MERLIN epochs reported here in which the 
source was resolved. The contour levels are indicated on the figure. 
Contours are at levels of (-2, -1.41, -1, 1, 1.41, 2, 2.828, 4, 5.657, 
8, 11.31, 16) x0.2 mjy / beam. The core and jet peak fluxes in 
this image are 1.24 0.06 and 1.18 0.06 mjy / beam respectively; 
the values listed in table 1 are for the entire of epoch Ml. 



at 2004 Feb 23 10:00. While available Knockin was included, 
although noisy, to improve the (u, v) plane coverage. 

The data reduction was done in the standard way for 
continuum mapping, outlined in the MERLIN user guide . 
We used 3C 286 as a flux density and polarisation position- 
angle calibrator, OQ 208 or 0552+398 as a point source 
calibrator. The approximately 400 mjy source 1951+355, 
which is unresolved to MERLIN at this frequency, was used 
as a phase reference and to calibrate the d-terms. The useful 
bandwidth was 14 MHz after editing the end channels. 

Amplitude and phase self-calibration was performed on 
the phase reference source (shifting it to the phase-centre) 
and interpolated to the target, Cygnus X-l. Observations 
taken on Cygnus X-l and 1951+355 were nagged at eleva- 
tions beneath 15° as phase and amplitude calibration trans- 
fer is generally poorer at low elevations. The imaging was 
performed using IMAGR within AIPS (Diamond 1995). A pre- 
liminary image showed that our coordinates used in the cor- 
relation for either the phase-reference or target source were 
slightly in error (~30 mas north-east of the actual source 
position for Cygnus X-l) and the phase centre was shifted 
with UVFIX within AIPS to the target position. As the proper 
motion of Cygnus X-l on timescales of months is very small 
compared to our resolution this shift was applied to all of 
the epochs. From these observations our best coordinates 
for Cygnus X-l are (J2000) RA 19 58 21.6747 Dec 35 12 
05.767; the peak core flux densities listed in table 1 were all 
measured within 4 mas of this position. 

All of the images were restored with a 50 mas circular 
beam. Table 1 summarizes the results from each image. All 
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Epoch 


Date (2004) 


Core peak 


Jet peak 






(mjy bm _1 ) 


(mjy bm ) 


Ml 


Jan 15/16 


1.2 + 0.07 


1.2 + 0.07 


M2 


Jan 26 


11.8 + 0.17 


< 0.51 


M3 


Feb 21 


9.9 + 0.20 


< 0.60 


M4 


Feb 22 


6.9 + 0.28 


< 0.84 


M5 


Feb 23 


7.8 + 0.30 


< 0.90 


M6 


Feb 25 


13.6 + 0.13 


< 0.39 



Table 1. Summary of the MERLIN observations of Cygnus X-l. 



epochs apart from epoch Ml, were found to be point sources 
(limits on the flux from jets are given in Table 1). No linear 
polarization was detected from the source at any time; the 
most stringent limit being P/I < 2.3% (3cr). 

In Fig 2 we present the image from epoch Ml, the only 
epoch at which a resolved jet was detected, and also the 
epoch with by far the weakest core emission. The source 
is resolved towards the north-west with an apparent bright 
component, possibly superposed on a smooth underlying jet. 
Applying the AIPS task JMFIT we find two approximately 
beam-sized Gaussian components of integrated intensities 
1.6 mjy (core) and 0.7 mjy (jet); see table 1. The position 
angle between the centres of these Gaussians is -25° ± 3°. 
The angular separation between core and jet is 70 + 5 mas 
(corresponding to 0.82 ± 0.06 light days at a distance of 2 
kpc). 

For comparison, epoch M6 is our image with the high- 
est dynamic range after the 120 mjy flare of 2004 Feb 20. 
We have used an identical beam size and contouring scheme 
as used for epoch Ml. There is no trace of extension in any 
position angle. To stress the point, had the extended 'jet' ob- 
served at epoch Ml been present in any of the other epochs, 
it would have been detected - its detection at epoch Ml is 
not just a consequence of the weaker core at this time. 

In an attempt to pin down the time of appearance of the 
extended radio emission, we have broken epoch Ml, which 
was a ~ 23 hr run from start to finish, into different seg- 
ments. A significant gap in the observations, between 2004 
Jan 15 21:00 - Jan 16 03:00 separates the data into two seg- 
ments, albeit of different length (6 and 15 hr respectively). 
The jet was only detected in the second of the two blocks of 
data. However, since MERLIN does not have enough base- 
lines to operate as a 'snapshot' type of instrument (unlike 
the VLA), the iw-coverage of the two blocks would have 
been very different. We therefore tested whether the jet ap- 
pearing in the second half of the data would have been de- 
tected by the «u-coverage of the first half. The clean com- 
ponents from the jet in the image of the second half were 
selected using the AIPS task CCEDT and added back in to the 
uv-Aa±a, of the first half using the task UVSUB. Imaging this 
modified tti>-data set, it was not possible to unambiguously 
detect the jet above the noise. A parallel approach was used 
in MIRIAD (Sault, Teuben & Wright 1995) with the task 
UVMODEL, and the same conclusions reached. From these 
data, we cannot therefore place a limit on the jet velocity 
based on epoch Ml alone. 
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Figure 3. Radio and X-ray monitoring around the time of the 
MERLIN observations. The top panel presents the Ryle Telescope 
15 GHz flux densities and the total (2-12 keV) X-ray count rate 
from the RXTE ASM. Arrows indicate the times of the MER- 
LIN observations; only the first observation (epoch Ml, see text) 
detected an extended jet. The dotted triangle in the upper panel 
indicates the time of the large radio flare plotted in Fig 4, which 
corresponded to a large increase in the X-ray flux. The lower panel 
indicates the hardness ratio from the ASM (5—12 keV/1.5— 3 keV), 
and the vertical lines indicate our best estimate of the range of 
launch times for the jet which brightened during epoch Ml. 
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Figure 4. A close up of the radio flare observed at 15 GHz with 
the RT on 2005 Feb 20, with 32-sec sampling. This is the highest 
radio flux density recorded from the source with the RT to date. 
MERLIN epochs M3-6 were performed between 1-5 days after 
this event, and yet none of them resolved a discrete ejection event, 
unlike epoch Ml. 



3 DISCUSSION 

3.1 Disc-jet coupling 

Fender, Belloni & Gallo (2004) have outlined a model in 
which discrete radio events are associated with transitions 
from 'hard' X-ray spectral states to 'softer' states at rela- 
tively high (^ 1% Eddington) luminosities. In the accretion 
flow this seems to be associated with a rapid inwards tran- 
sition by the optically thick, geometrically thin, accretion 
disc. In the outflow it appears to correspond to a transi- 
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Figure 5. A hardness intensity diagram for Cyg X-l, based upon 
the RXTE ASM one-day averages. The epochs of the MERLIN 
observations, previous VLBA observations which revealed a com- 
pact (^20 mas) steady jet, and the large radio flare observed by 
the RT (Fig 4), are indicated. The observation of the jet is the 
only one to the left of the peak which we might tentatively asso- 
ciate with the 'jet line' in this source. However, this interpretation 
does not explain the nature of the 2004 Feb 20 flare (Fig 4) , which 
appears to have occurred to the right of the line, yet we would 
naturally interpret as being associated with an ejection / shock 
event. 'HIMS' = 'Hard Intermediate State'; 'SIMS' = 'Soft In- 
termediate State'; 'SPL' = 'Steep Power Law' - see Fender et al. 
(2004); Homan & Belloni (2005); McClintock & Remillard (2005); 
Belloni (2005). 



tion between a state which produces a steady jet and one 
which produces no jet, via a brief phase in which the jet 
Lorentz factor increases, resulting in internal shocks in the 
flow. It may be these internal shocks which we identify as 
discrete components in radio maps of X-ray transients (see 
also Kaiser et al. 2000). 

The behaviour of Cyg X-l appears to follow a consistent 
pattern, although with some anomalies. Figure 3 presents a 
close-up of the data around the time of the one resolved 
event, 2004 Jan 15/16. Precisely around the time of the 
MERLIN observation the radio flux dropped to very low 
levels, towards the end of a period of X-ray spectral soften- 
ing. The radio flux stayed at these very low (< 2 mjy at 5 
GHz and undetectable at 15 GHz) levels for 12-18 hr and 
then recovered rapidly recovered to previous levels 6-12 hr 
before the X-ray spectrum returned to a hard state. Unfor- 
tunately the RXTE ASM sampling was rather sparse around 
this time and this hampers our precise interpretation of the 
sequence of events. 

Fig 3 also shows that the radio flare of 2004 Feb 20 
(Fig 4) occurred around the time of a large step-like in- 
crease of the integrated X-ray flux, while the source was 
in a 'hard intermediate' state (Fender et al. 2004; Homan 
& Belloni 2005). Fig 5 (see also Fig 6) sheds more light 
on this in the form of a hardness-intensity diagram (HID) 
for Cyg X-l. There is a clear peak in the luminosity dis- 
tribution of the source, corresponding to a hardness ratio 
HR2 (counts in 5-12 keV band / counts in 1.5-3 keV band) 
around 0.3-0.4. In the canonical hard state 1.0 <HR2< 1.5 
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Figure 6. Motion of Cyg X-l in the HID for four days preceding 
and two days following the jet which was imaged at epoch Ml. 



is more typical. It is interesting to note that epoch Ml - the 
only epoch at which a jet was imaged - is the only epoch 
to the left of (softer than) the peak. All of the other MER- 
LIN epochs correspond to 'hard intermediate' states (HIMS; 
Belloni 2005) with 0.5 <HR2< 1.0. We therefore tentatively 
associate HR2~ 0.4 with the 'jet line' (Fender et al. 2004) 
for Cyg X-l. This is consistent with the picture of Fender et 
al. (2004) and Corbel et al. (2004) in which the steady jet 
(core emission) persists during softening from the canonical 
hard state, until some critical point (represented spectrally 
by the 'jet line') is reached. Also marked in Fig 5 are the 
state of Cyg X-l at epochs when a steady, compact (^20 
mas) jet was imaged with the VLBA (Stirling et al. 2001; R. 
Spencer private communication); these are all in the canon- 
ical low/hard state. 

Fig 6 shows the evolution of daily- averaged counts for 
Cyg X-l in the HID for the four days preceding and two 
days following the imaging of the jet on 2004 Jan 15/16. 
Based on this more detailed analysis, we find that the jet 
line is likely to lie close to HR~0.3; this is indicated on the 
figure. We note that Cyg X-l exhibits X-ray temporal and 
spectral variability on much shorter timescales than those 
utilised here from the RXTE ASM, and we may be miss- 
ing key events by using daily averages. It is also clear that 
Cyg X-l does not behave like a 'typical' transient once it 
has crossed the jet line, in progressing towards the canon- 
ical soft / thermal state, but instead returns to the hard 
IS state within two days. The only other source to exhibit 
such behaviour is GRS 1915+105, although it is likely that 
repeatedly radio-bright objects such as Cygnus X-3 may be- 
have similarly. 

It is also interesting to speculate about the nature of 
the core emission at epoch Ml. It has been well established 
for some years now (e.g. Fender et al. 1999) that the radio 
emission from black hole X-ray binaries is suppressed in soft 
X-ray states compared to hard states. An unambiguous de- 
tection of core (i.e. currently generated flow) radio emission 
from a source in a steady soft state has not yet been re- 
ported. If the 'core' radio emission we detect at epoch Ml is 
indeed the supressed level, then the quenching in Cyg X-l 
is only a factor of ~ 10, compared to > 30 and > 50 in the 
black hole X-ray binaries GX 339-4 and XTE J1550-564 re- 
spectively (Fender et al. 1999; Corbel et al. 2001). However, 



there is no way of knowing if this weak core is completely 
or only partially quenched given, again, the known rapid 
variability of Cyg X-l. What we observe as the core could 
in fact be emission from a receding jet which in many cases 
will have a very small proper motion (much less than that of 
the approaching jet). In this scenario the jet emission would 
have to have some flux evolution including a local peak; for 
a symmetric monotonically decaying event the receding jet 
cannot appear brighter than the approaching. 



3.2 Jet launch and velocity 

As discussed above, it is very difficult to pin down the mo- 
ment at which point the matter / energy associated with the 
resolved jet was launched. In the model of Fender, Belloni 
& Gallo (2004) the launch moment would correspond to the 
end of a phase of X-ray spectral softening and would be im- 
mediately followed by a quenching / suppression of the core 
radio emission. However, given the limited quality of data 
from the RXTE ASM, even for a very bright source such 
as Cyg X-l, we can only say that < At < 15 where At 
was the time since launch in days. The resulting associated 
proper motion and velocity would be fi — 70/ At mas d _1 
and Wapp = 0.8c/At. The more detailed inspection of the 
HIDs outlined above suggests - no more - that the jet was 
launched at most ~ 4 days before epoch Ml. This interpre- 
tation is supported circumstantially by the very rapid decay 
of the large flare event on 2005 Feb 20 (Fig 4). Further- 
more, we estimate that proper motions larger than ~ 100 
mas d _1 would probably have caused smearing in the Ml 
image which was not observed. 

For 0.5 < At < 4 days, 140 < fi < 17.5 mas cT 1 , 
1.6c < i> app < 0.2c. The Cyg X-l system has been exten- 
sively modelled based on optical and UV data, from which 
a general consensus has emerged that the orbital inclination 
angle of the system is ~ 30° (e.g. Gies & Bolton 1986b). 
If the jet shares the same inclination angle (and this is by 
no means certain - Maccarone 2002) then the true velocity 
of the jet would be ^ 0.3c (given the accumulation of un- 
certainties, we cannot place an upper limit on the Lorentz 
factor of the jet - see Fender 2003; Miller- Jones, Fender & 
Nakar 2005 for further discussion). 



4 CONCLUSIONS 

We report the first observation of a transient, extended (70 
mas) radio jet from the classical black hole binary Cygnus 
X-l. The jet was observed in only one of six MERLIN obser- 
vations performed while the system was undergoing frequent 
X-ray state changes. The one epoch in which the jet was im- 
aged was that with the softest 1.5-12 keV spectrum; at that 
moment the core radio emission was detectable but a factor 
~ 10 weaker than normal. Estimates of the launch epoch 
of the jet indicate it probably has a bulk velocity v > 0.3c, 
with essentially no formal upper limit. 

At the other five epochs the source X-ray spectrum was 
harder, although still not in the canonical hard state, and 
displayed more powerful core emission and no extended jet. 
Three of these observations were made within days of the 
brightest radio flare observed from Cyg X-l (~ 120 mjy 
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at 15 GHz), and yet no extended radio emission was de- 
tected. It is not clear if this may be related to the jet ex- 
panding into a rarefied bubble (Gallo et al. 2005), or if for 
some reason the flare was, unusually, not associated with 
a jet. Throughout the interpretation of the X-ray:radio be- 
haviour, we have to bear in mind that Cyg X-l varies on 
much shorter timescales than the RXTE ASM sampling, 
and we may have missed some key event. Nevertheless, the 
overall pattern of behaviour seems to be broadly consistent 
with the 'unified' model for black hole X-ray binary jets put 
forward by Fender, Belloni & Gallo (2004) although there 
are some clear anomalies. 
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